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Abstract. By means of the multifractal analysis (MFA), the expressions of the probability density 
functions (PDFs) are unified in a compact analytical formula which is valid for various quantities in 
turbulence. It is shown that the formula can explain precisely the experimentally observed PDFs 
both on log and linear scales. The PDF consists of two parts, i.e., the tail part and the center 
part. The structure of the tail part of the PDFs, determined mostly by the intermittency exponent, 
represents the intermittent large deviations that is a manifestation of the multifractal distribution of 
singularities in physical space due to the scale invariance of the Navier-Stokes equation for large 
Reynolds number. On the other hand, the structure of the center part represents small deviations 
violating the scale invariance due to thermal fluctuations and/or observation error. 



INTRODUCTION 

In this paper, we derive the unified formula for various probability density functions 
(PDFs) in fully developed turbulence by means of the multifractal analysis (MFA) HHH 
SIII&IEBIESEiIEIEIIK and analyze the PDFs observed in two experiments, i.e., the 
PDFs of velocity fluctuations, of velocity derivatives and of fluid particle accelerations at 
R\ = 380 that was extracted by Gotoh et al. from the DNS of the size 1024 3 JH], and the 



PDF of fluid particle accelerations at R\ = 690 obtained in the Lagrangian measurement 
of particle accelerations that was realized by Bodenschatz and co-workers ful [l5l Hill 
by raising dramatically the spatial and temporal measurement resolutions with the help 
of the silicon strip detectors. The MFA of turbulence is a unified self-consistent approach 
for the systems with large deviations, which has been constructed based on the Tsallis- 
type distribution function H17L ll 811 that provides an extremum of the extensive Reny 
lfl9ll or the non-extensive Tsallis entropy tl 1VL fTsL I20I1 under appropriate constraints. The 
analysis rests on the scale invariance of the Navier-Stokes equation for high Reynolds 
number, and on the assumptions that the singularities due to the invariance distribute 
themselves multifractally in physical space. The MFA is a generalization of the log- 
normal model lEILl22Ll23ll . It has been shown |@] that the MFA derives the log-normal 
model when one starts with the Boltzmann-Gibbs entropy. 

For high Reynolds number Re ^> 1, or for the situation where effects of the kinematic 
viscosity v can be neglected compared with those of the turbulent viscosity, the Navier- 
Stokes equation, du/dt + (u • V)u = —V (p/p) + ^V 2 u, of an incompressible fluid is 



invariant under the scale transformation Il24ll25ll2al r — > Ar, u — > A a / 3 u, £ — > A 1 a ^t and 



(p/p) ~~ * A 2a//3 (p/p) where the exponent a is an arbitrary real quantity. The quantities 
p and p represent, respectively, mass density and pressure. The Reynolds number Re 
of the system is g iven by Re = 5u^£- m /v = (£- m /vi) 4 ^ w i m the Kolmogorov scale 
77 = (z/ 3 /e) 1/ ' 4 12711 where e is the energy input rate at the input scale £- in . Here, we 
introduced 5u- in = \u(» + £- m ) — u(»)\ with the definition of the velocity fluctuation 
(difference) 5u n = \u(» + £ n ) — it(«) | where u is a component of velocity field u, and £ n 
is a distance between two points. The pressure (divided by the mass density) difference 
Sp n = \p/ p(» + £ n ) —pi p(«) I between two points separated by the distance £ n is another 
important observable quantity. We are measuring distance by the discrete units £ n = S n £o 
with 5 n = 2~ n [n = 0,1,2, •••). The non-negative integer n can be interpreted as the 
multifractal depth. However, we will treat it as positive real number in the analysis of 
experiments. The multifractal depth n is related to the number of steps within the energy 
cascade model. 1 At each step of the cascade, say at the nth step, eddies break up into 
two pieces producing an energy cascade with the energy-transfer rate e„ that represents 
the rate of transfer of energy per unit mass from eddies of size £ n to those of size £ n+ \. 



SINGULARITIES AND SCALING EXPONENTS 

Let us consider the quantity 5x n = \x(* + £ n ) — x(»)\ having the scaling property 
\x n \ = \Sx n /5xo\ = 5% a / 3 . Its spatial derivative is defined by \x'\ = lini£ n ^o 5x n /£ n oc 
Hindoo l^/ 3 " 1 which becomes singular for a < 3/0. The values of exponent a specify 
the degree of singularity. We see that the scale invariance provides us with 5u n /Suo = 
and 5p n /5p = (£ n /£o) 2a ^ 3 giving, respectively, = 1 for the velocity fluctuation 
and 4> = 2 for the pressure fluctuation. The velocity derivative and the fluid particle ac- 
celeration may be estimated, respectively, by \u'\ = Hindoo and by |a| = Hindoo a n 
where we introduced the nth velocity derivative u' n = 5u n /£ n and the nth fluid particle 
acceleration a n = Sp n /£ n corresponding to the characteristic length £ n . Note that the ac- 
celeration a of a fluid particle is given by the substantive time derivative of the velocity: 
a = du/dt + (u ■ V)u. We see that the velocity derivative and the fluid particle accelera- 
tion become singular for a < 3 and a < 1.5, respectively, i.e., \u'\ oc lim^^o^/ 3 ^ 1 — > 
00 and |a| oc lini£ n ^o4 l 2a ^ 3 ' ) 1 ~" * 00 • We also see that the energy dissipation rate be- 
comes singular in the limit n — ► 00 for a < 1, i.e., Hindoo e„/e = lim„^ oo (£ n /£ ) a_1 — > 
00 giving = 3. 

The MFA rests on the multifractal distribution of singularities that is a manifes- 
tation of the scale invariance of the Navier-Stokes equation for large Re as men- 
tioned above. The probability P^(a)da to find, at a point in physical space, a sin- 
gularity labeled by an exponent in the range a ~ a + da is given by H, 0, @] 



1 The definition of the number of steps n within the energy cascade model is given by n — — log 2 (r/£ in ) 
for the eddies whose diameter is equal to r. By putting r = £ n , this gives us the relation between n and n 
in the form 

n = n-log 2 (£ /e- m ). (1) 




FIGURE 1. Comparison of the present scaling exponents £ m for ji = 0.238 (solid curve) with the 
experimental results plotted by circles at R\ = 110 (Re — 32000) [30], and with other theories with 
the same value of /i. K41 is given by the dotted line, /3-model by the dashed line, p-model by the dotted 
dashed line, log-Poisson model by the short dashed curve, and log-normal by the two dotted dashed curve. 



p( n \a) = [1 — (a — a )V(^ a ) 2 ] /^i™^ w ith an appropriate partition function 

and {Aaf = 2X/[(1 - q) In 2] . This is consistent with the relation JHH ?W (a) oc 
8]^^°^ that reveals how densely each singularity, labeled by a, fills physical space. 
In the present model, the multifractal spectrum f(a) is given by 0,13,0.13] /(a) = 
1 + (1 — g) _1 log 2 [l — (a — a ) 2 /(Aa) 2 ]. The range of a is a min < a < a max with 
ttmin = olq — Aa, a max = «o + Aa. The distribution function p( n \a) is determined by 
taking an extremum of generalized entropies, i.e., the extensive Renyi entropy or the 
non-extensive Tsallis entropy, under the condition that the information one has for the 
system is only the value of the intermittency exponent. In spite of the different charac- 
teristics of the entropies, the distribution functions (a) giving their extremum have 
the common structure. 2 

The dependence of the parameters «o, X and q on the intermittency exponent fi is 
determined, self-consistently, with the help of the three independent equations, i.e., the 
energy conservation: (e n ) = e, the definition of the intermittency exponent fi: (e 2 ) = 
e 2 5~ M , and the scaling relation 3 : 1/(1 — q) — 1/ct- — l/a + with a± satisfying f(a±) = 0. 
The average (• • ■) is taken with P^ n \a). 

The scaling exponents Cm of the mth order velocity fluctuations, defined by (\u n \ m ) = 
( 5™ a/3 } oc , are given in the analytical form H 5 H 

( m = a m/3 - 2Xm 2 /[9(l + Cj£ )] - [1 - log 2 (l + C 1 ^)]/^ - q) (2) 

with Cq = 1 + 2g 2 (l — q)Xln2. The formula © is independent of n, that is a manifes- 
tation of the scale invariance. 



Within the present formulation, the decision cannot be pronounced which of the entropies is underlying 
the system of turbulence. 

3 The scaling relation is a generalization of the one derived first in i28l l29il to the case where the 
multifractal spectrum has negative values. 



The derived scaling exponents © are shown in Fig. [U by the solid curve for the case 
[i = 0.238, and are compared with experimental data [ 30] and with the curves given by 
other theories, i.e^pi log-normal JULIHIIl, /3-model iU, p-model TH Ej] 
and log-Poisson SH- 



VARIOUS PROBABILITY DENSITY FUNCTIONS 

It has been shown that the probability Tl^\(x n )dx n to find a physical quantity x n in the 
range x n ~ x n + dx n is given in the form 

U^\x n )dx n = U^(x n )dx n + AU^\x n )dx n (3) 

with the normalization ff^ dx n TvJ/ (x n ) = 1. The first term represents the contribution 
by the singular part of the quantity x n stemmed from the multifractal distribution of its 
singularities in physical space. This is given by II^ s (|:r n |)<i;E n oc P^(a)da with the 

transformation of the variables, |sc n | = 5^ Q//3 . Whereas the second term AH^\x n )dx n 
represents the contribution from the dissipative term in the Navier-Stokes equation, 
and/or the one from the errors in measurements. The dissipative term has been dis- 
carded in the above investigation for the distribution of singularities since it violates the 
invariance under the scale transformation. The contribution of the second term provides 
a correction to the first one. Note that each term in © is a multiple of two probability 
functions, i.e., the one to determine the portion of the contribution among the above 
mentioned two independent origins, and the other to find x n in the range x n ~ x n + dx n . 
Note also that the values of x n originated in the singularity are rather large representing 
intermittent large deviations, and that those contributing to the correction terms are of 
the order of or smaller than its standard deviation. 

The mth moment ofthe variable \x n \ is given by {{\x n \ m ))^ = f^dx n \x n \ m Ily (x n )= 

2 7 <2, + (1 - 2 7 $)a 0m &~ where 2 7 £ = J^dx^x^AU^M, and a 3 , = 

{2/^(1 + v / aJ)]} 1 / 2 . 

We now derive the PDF, fl^ (£ n ), defined by the relation flf ] (f n )df„ = uf } (x n )dx n 

with the variable £ n = x n j '((l^nl 2 ))^ 2 normalized by the standard deviation ((x n ))^ 2 . 
This PDF is to be compared with the observed PDFs. The variable is related with a 
by |£n| = Zjt /3 ~ W2 with £„ = [2 7 g5n C20 + (1 - 2 7 $)a 20 ]- 1 /2. ft is reasonable to 
imagine that the origin of intermittent rare events is attributed to the first singular term 
in ©, and that the contribution from the second term is negligible. We then have for the 
tail part, i.e., C < IU < C ax > 

ni n) (£nR., 



fi-(n) £n 1 

\U 



l-q (31n|e„/e 
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20 2 X|ln5 n | 
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FIGURE 2. Analyses of the PDFs of velocity fluctuations (closed circles) and of velocity derivatives 
(open circles) measured in the DNS by Gotoh et al. at R\ = 380 by the present theoretical PDFs IT") (£ n ) 
for velocity fluctuations (solid lines) and for velocity derivatives (dashed line) are plotted on (a) log and 
(b) linear scales. The DNS data points are symmetrized by taking averages of the left and the right 
hand sides data. The measuring distances, r/77 = i n /rj, for the PDF of velocity fluctuations are, from 
the second top to bottom: 2.38, 4.76, 9.52, 19.0, 38.1, 76.2, 152, 305, 609, 1220. For the theoretical PDFs 
of velocity fluctuations, fi = 0.240 (q = 0.391), from the second top to bottom: (n, n, q') = (20.7, 14.6, 
1.60), (19.2, 13.1, 1.60), (16.2, 10.1, 1.58), (13.6, 7.54, 1.50), (11.5, 5.44, 1.45), (9.80,3.74, 1.40), (9.00, 
2.94, 1.35), (7.90, 1.84, 1.30), (7.00, 0.94, 1.25), (6.10, 0.04, 1.20), £* = 1.10 - 1.43 (a* = 1.07), and 
£f ax = 204 ~ 6.63. For the theoretical PDF of velocity derivatives, (n, n, q') = (22.4, 16.3, 1.55), 
I* = 1.06 (a* = 1.07), and £™ ax = 302. For better visibility, each PDF is shifted by -1 unit along the 
vertical axis. 



with n 



in) 



3(1 - 2 1 ! ) n) )/(2 ( f ) i n J27cX\ln5 n \), £ n . 



^^ Qmin/3_C2 * /2 . On the other hand, for the center part, the contribution to the PDF 
comes, mainly, from thermal fluctuations or measurement error. It may be described by 
the Tsallis distribution function with respect to the variable £ n itself, i.e., |£ n | < £*, 



ni n) (£nR, 




1/(1-9') 



<% n . (5) 



This specific form of the Tsallis distribution function is determined by the condition that 
the two PDFs © and © should have the same value and the same slope at £* which is 

defined by C = Ut* /3 ~ C ^ /2 with a* being the smaller solution of C20/2 - (pa/3 + 1 - 

f(a) = 0. It is the point at which fl^ (£*) has the least n-dependence for large n. 
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FIGURE 3. Comparison between the PDF of fluid particle accelerations measured in the DNS by 
Gotoh et al. at R\ = 380 and the present theoretical PDF A^(u n ) are plotted on (a) log and (b) 
linear scales. Closed circles are the DNS data points both on the left and right hand sides of the PDF. 
Solid lines represent the curves given by the present theory with /i = 0.240 (q — 0.391), (n, n, q') — 
(17.5, 11.4, 1.70), lj* = 0.622 (a* = 1.01), and < nax = 2530. 



With the help of the second equality in ©, we obtain AIT^ (x n ), and have the formula 
to evaluate 7 £ in the form 2 7 ^ = (tf - ) / (l + K<$ - L$) where 
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with = ^min/C = <5^ Qmax - Q *)/ 3 , 2* = infl/H = <J*(°o-«')/3 '. We see that the tail part 
of the PDF, given by ©, is mostly determined by the intermittency exponent \i and the 
multifractal depth n which gives a length scale i n . On the other hand, the center part of 
the PDF, ©, is mainly controlled by q'. 

The PDFs both for velocity fluctuations and for velocity derivatives are given 
by the common formula f[( n \£ n ) = n^ 1 (^ n ) in their normalized variables 
£n = b~Un/ (((dun) 2 }} 1 / 2 . On the other hand, the PDFs both for pressure differences and 
for fluid particle accelerations are given by the common formula A^(uj n ) = f[^ 2 (a; n ) 

in their normalized variables u n = 5p n / ' (((Spn) 2 )) 1 ^ 2 ■ The PDF for energy dissipation 
rates is given with = 3. 

The PDFs extracted by Gotoh et al. from their DNS data JH at R x = 380 are 
shown, on log and linear scales, in Fig.|2]both for velocity fluctuations and for velocity 
derivatives, and in Fig.^for fluid particle accelerations. We found the value \i = 0.240 
by analyzing the measured scaling exponents ( m of velocity structure function with the 
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FIGURE 4. Dependence of q' on the distance r/rj extracted from the analyses of the PDFs for velocity 
fluctuations (closed circles), for velocity derivatives (open circle) and for fluid particle accelerations (open 
triangle). The line represents q' — — 0.051og 2 (r/ry) + 1.71. 

formula ©, which gives the values q = 0.391, a = 1.14 and X = 0.285. Through the 
analyses of the PDFs for velocity fluctuations in Fig. O we extracted the formula for the 
dependence of n on r/77: 0,0] 

n = -0.989 xlog 2 r/r/+ 16.1 (for4<r), (8) 
n = -2 AO x\og 2 r/r] + 24.0 (forr<4). (9) 

This shows that the inertial range is divided into two scaling regions separated by the 
characteristic length £ c /r) = 48.7 which is close to the Taylor microscale A/77 = 38.3 of 
the system. The equation dU) is consistent with the picture of the energy cascade model 
in which each eddy breaks up into 2 pieces at every cascade steps, whereas © indicates 
that, for r < £ c , each eddy breaks up, effectively, into 1.33 ~ 4/3 JSJ pieces at every 
cascade steps. This fact may be attributed to a manifestation of structural difference of 
eddies, which can be checked by visualizing DNS eddies. Actually, one observes that 
DNS eddies with larger diameters than Taylor microscale A have rather round shapes, 
whereas eddies with smaller diameters compared with A have rather stretched shapes 
ll35ll . The energy input scale for this DNS is estimated as the longest scale available in the 
lattice with cyclic boundary condition, i.e., im/r] = 71-/77 ~ 1220 with 77 ps 0.258 x 10~ 2 
lfl3ll which gives the number of steps n within the energy cascade model through the 
formula © with i /r] w 81300 determined by (El). 

For the analysis of the PDF for velocity derivatives in Fig. |2l we chose the value 
(n, n, q') = (22 A, 16.3, 1.55). The length corresponding to n is calculated by © to 
give r/77 = 1.61, which may provide us with an estimate for the effective shortest length 
in processing the DNS data to extract velocity derivatives. Note that it is about the same 
order of the mesh size Ar/r) = 2tt/(1024 x 77) w 2.38 [13] of the DNS lattice. 

For the PDFs for fluid particle accelerations in Fig. El we have (n, n, q') = 
(17.5, 11.4, 1.70). Substitution of this value into ® gives the corresponding character- 
istic length r/77 = 7.91 IIHD - This may be the effective minimum resolution in cooking 
the DNS data to distill accelerations. 
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FIGURE 5. Comparison between the experimentally measured PDF of fluid particle accelerations by 
Bodenschatz et al. at R\ = 690 (Re = 31 400) and the present theoretical PDF A (n) (u n ) are plotted on (a) 
log and (b) linear scales. Open squares are the experimental data points on the left hand side of the PDF, 
whereas open circles are those on the right hand side. Solid lines represent the curves given by the present 
theory with fi = 0.240 (q = 0.391), (n, q') = (17.1, 1.45), w* = 0.605 (a* = 1.01), and w™ ax = 2040. 

In Fig. we plotted the dependence of q' on the distance r/t] extracted from the 
analyses of the DNS data, i.e., the closed circles are extracted from the PDF of velocity 
fluctuations, the open circle is from the PDF of velocity derivatives and the open triangle 
from the PDF of fluid particle accelerations. The line represents q' = — 0.051og 2 (r/^) + 
1.71. The points for r/77 > 20 (closed circles) and for the accelerations (open triangle) 
are quite sensitive and easy to be determined. Other points are insensitive and have a 
rather wide range in deciding the values q' . 

The PDFs for fluid particle accelerations measured by Bodenschatz et al. at R\ = 690 
iflfjl are given in Fig.|5J We determined the value n = 17.1 by substituting the reported 
value £ = £ in =7.1 cm and the spatial resolution 0.5 fim of the measurement for i n into 
its definition, n = log 2 (^o/^n)- The values /i = 0.240 and q' = 1.45 are extracted by the 
analysis of the experimental PDF with the derived theoretical formula 111 ill . Then, we 
have the values of parameters: q = 0.391, «o = 1-14 and X = 0.285. The flatness of the 
PDF turns out to be JHI] = ((a*)) / '((a^)) 2 = ((£)) = 56.9 which is compatible with 
the value of the flatness ~ 55 ± 4 reported in iflfill . 



DISCUSSIONS AND PROSPECTS 

In this paper, the various experimental PDFs in turbulence are analyzed precisely with 
the formulae © and © of the corresponding PDFs derived by the MFA. It is revealed 
that there are two distinct mechanisms underlying the dynamics of turbulence. One 
contributes to the tail part of PDFs and the other to the center part. The structure of the 
tail part of the PDFs is determined by the global structure representing the intermittent 
large deviations that is manifestations of the multifractal distribution of singularities 



in physical space due to the scale invariance of the Navier-Stokes equation for large 
Reynolds number. The specific form of the tail part comes from the assumption that the 
probability to find a singularity exponent a within the range a ~ a + da at a point 
in physical space is given by the Tsallis-type distribution function with the Tsallis 
parameter q. The relation between a and an observing variable is given by the scale 
transformation. On the other hand, the structure of the center part represents small 
deviations violating the scale invariance due to thermal fluctuations and/or observation 
error. The center part is assumed to be given by the Tsallis-type distribution function 
with the Tsallis parameter q' for the observing variable itself. The value of q' may be 
determined by a local structure of the system, e.g., the dynamics of a vortex, the mutual 
interaction between vortices and so on, and depends on the distance of two measuring 
points in contrast to q. The latter parameter q does not depend on the distance, and 
is determined once the value of the intermittency exponent is given. It is one of the 
attractive future problems to derive two different dynamics taking care of the tail part 
controlled by q and the center part by q', and will be reported in the near future. 

The success of the MFA in analyzing turbulence in high accuracy may provide us with 
a good tool to see what is the origin of the singularities and why their distribution is mul- 
tifractal. In order to investigate them, the vortex tangle is one of the attractive candidates 
as Feynman proposed [36], since the vorticity in superfluid 4 He and 3 He is quantized and 
the normal component within the sense of the two fluid model can be negligible at very 
low temperature. If the singularity originates from the core of vortex, the multifractality 
of turbulence in normal fluid can be related to various values of vorticities in the fluid. In 
this case, the vortex tangle may be uni-fractal, and does not exhibit intermittency. If the 
singularity originates from the reconnection of vortices, the multifractality of turbulence 
in normal fluid is related to the distribution of reconnection points in the fluid. Then, the 
vortex tangle may be also multifractal, and does exhibit intermittency. A temperature- 
independent vortex decay mechanism below T ~ 70 mK has been observed in superfluid 
4 He S, and the Kolmogorov spectrum (K41) is extracted from the simulation within 
the vortex filament model for non-frictional suferfluid 4 He ll3^1 . In tangle, the quantized 
vortex crosses the stream lines of superfluid velocity field, which may result in the decay 
mechanism at low temperature. We expect that, through the analysis of the local dynam- 
ics controlled by the Tsallis parameter q', we can extract some information about the 
mutual friction between the superfluid and normal components. 

A circular vortex lattice is observed in the simulation of fast rotating Bose-Einstein 
condensate confined in a 2-dimensional quadratic-plus-quartic potential [39]. A possi- 
bility of generation of tangle phase in this system is one of the attractive problems. The 
situation may be quite similar to the one in the generation of the Taylor-Couette turbu- 
lence in normal fluid. 

Let us close this paper by mentioning that the vortex tangle can be an important 
stationary phase of quLSI (quantum LSI) consisting of, for example, a huge number of 
superconducting loops of flux qubit [40]. The direction of current in a loop may change 
so quickly under an operation of quantum computer. The congeries of the flux qubits can 
produce a tangle phase of flux quantums. It may be important to see if the tangle phase 
benefits quantum entangled states or not. 
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